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The superconducting state of a thin superconducting disk 
with a hole is studied within the non-linear Ginzburg-Landau 
theory in which the demagnetization effect is accurately taken 
into account. We find that the flux through the hole is not 
quantized, the superconducting state is stabilized with in- 
creasing size of the hole for fixed radius of the disk, and a 
transition to a multi- vortex state is found if the disk is sufh- 
ciently large. Breaking the circular summetry through a non 
central location of the hole in the disk favours the multi- vortex 
state. 
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I. INTRODUCTION 

Thanks to recent progress in microfabrication and 
measurement techniques, it is possible to study the prop- 
erties of superconducting samples with sizes comparable 
to the penetration depth A and the coherence length ^. 
The properties of these mesoscopic systems are consid- 
erably influenced by confinement effects. Therefore, the 
vortex state will depend on the size and the geometry of 
the sample. 

In the present paper, we study the properties and the 
vortex states of superconducting thin disks with a hole. 
In the past, two limiting cases were studied: the thin- 
wire loop and the disk without a hole. In 1962 Little and 
Parks studied a thin-wire loop in an axial magnetic field 
PI . The Tc — H phase diagram showed a periodic compo- 
nent. Each time a flux quantum (J}q = hc/2e penetrates 
the system, Tc{H) exhibits an oscillation. Berger and Ru- 
binstein |2| studied mesoscopic superconducting loops us- 
ing the non-linear Ginzburg-Landau (GL) theory. They 
assumed that the induced magnetic field can be neglected 
for samples with sufficiently small thickness. In the limit 
of thin loops, the transition between states with differ- 
ent angular momentum L (also called vorticity) occurs 
when the enclosed flux (j) equals {L + 1/2) ^o [||. The su- 
perconducting disk was studied by Schweigert et al 0-R] 
(see also Ref. [||) by solving the two GL equations self- 
consistently. Although the GL equations were derived to 
describe superconductivity near the critical point, this 
theory turns out to be valid over a much broader range 
of magnetic field and temperature Mpl . They found that 
the finite thickness of the disk infiuences the magnetic 
field profile, i.e. the magnetic pressure, and this changes 
the size of the Meissner effect which is different from the 
well-studied cylinder geometries pO|]. The reverse prob- 



lem, i.e. the anti-dot, was studied by Bezryadin et al 
|ll| . They obtained a phase diagram of a thin super- 
conducting film with a circular hole in an axial magnetic 
field by solving numerically the nonlinear GL equations 
in the limit of a thin film. Here we generalize the re- 
sults of Ref. pI to a thin circular superconducting disk 
containing a circular hole. 

The intermediate case of finite width loops was stud- 
ied previously by Bardeen jl^l within the London theory. 
He showed that in tubes of very small diameter and with 
wall thickness of the order of the penetration depth the 
fiux through the tube is quantized in units of z^0o where 
V < 1 depends on the dimensions of the system. Arutun- 
yan and Zharkov (iqjljl found that the fiux through the 
effective area tt {p*) is equal to r7i(/)o, where the effective 
radius p* is approximately equal to the geometric mean 
square of the inner radius Ri and the outer radius Ro of 

1 /2 

the cylinder; i.e. p* = (RiRo) . Recently, Fomin et al 
||l5| studied square loops with leads attached to it and 
found inhomogcncous Cooper pair distributions in the 
loop with enhancements near the corners of the square 
loop. Bruyndoncx et al n& investigated infinitely thin 
loops of finite width. In this case, the magnetic field 
induced by the supercurrents can be neglected and the 
total magnetic field equals the external applied magnetic 
field. Furthermore, they used the linearized GL equa- 
tion which is only valid near the superconductor/normal 
boundary where the density of the superconducting con- 
densate IV' I is small. Only the giant vortex state with 
a definite angular momentum L was studied and they 
concentrated on the two- (2D) to three-dimensional (3D) 
crossover. Berger and Rubinstein [|l^ also studied in- 
finitely thin loops of finite width using the non-linear GL 
theory and they also neglected the induced field. 

It is well known that for type-II (k > l/V^) supercon- 
ductors the triangular Abrikosov vortex lattice is ener- 
getically favored in the range Hd < H < Hc2- Since the 
effective London penetration depth A = A^/d increases 
considerably in thin samples and for d <^ X one would ex- 
pect the appearance of the Abrikosov multi- vortex state 
even in thin rings made with a material with k < l/v2. 
Similar as for the case of a thin disk , we expect that 
the structure of a finite number of vortices should differ 
from a simple triangular arrangement in the case of a thin 
ring and we expect that they will resemble the configura- 
tions found for Coulomb clusters which are confined into 
a ring ^. 

In the present paper we consider circular fiat disks of 
non-zero width with a circular hole in it, which is not 



necessary in the center of the disk. The superconducting 
properties are also studied deep inside the superconduct- 
ing state where: i) nonUnear effects are important, i.e. 
1^1 is not necessarily small, and the nonlinear GL equa- 
tions have to be solved, ii) the total magnetic field is not 
homogeneous, i.e. it is spatial varying due to the Meiss- 
ner effect and the flux quantization condition which may 
enhance or diminish the magnetic field through the hole 
as compared to the applied magnetic field, and iii) due 
to nonlinear effects the circular symmetric giant vortex 
states are not necessarily the lowest energy states and the 
magnetic field can penetrate the superconductor through 
single vortices creating a multi-vortex state. 

The paper is organized as follows: fn Section II we 
present the theoretical model. In Section III we consider 
a small superconducting disk with a hole in the center. 
In this case we find that only the giant vortex state ap- 
pears. We study the influence of the radius of the hole 
on the superconducting state. For such a small system 
the relation between the local magnetic field, the current 
density and the Cooper pair density is investigated and 
the quantization of the flux through the hole is inves- 
tigated. Next, in Section IV, we consider the case of a 
larger superconducting disk with a hole in the center. For 
increasing magnetic field, we find re-entrant behaviour; 
i.e. transition from the giant vortex state to the multi- 
vortex state and back to the giant vortex state before 
superconductivity is destroyed. Finally, in Section V, we 
investigate the influence of the position of the hole on the 
vortex configuration; i.e. what happens if we break the 
axial symmetry? Our results are summarized in Section 
VI. 



II. THEORETICAL FORMALISM 

In the present paper, we consider superconducting 
disks with radius Ro and thickness d with a hole inside 
with radius i?, which is placed a distance a away from 
the center of the disk (Fig. p. These superconducting 
'rings' are immersed in an insulating medium with a per- 
pendicular uniform magnetic field Hq . To solve this prob- 
lem we follow the numerical approach of Schwcigert and 
Peeters [Q. For thin disks {d <$^ ^, A) they found that 
it is allowed to average the GL equations over the disk 
thickness. Using dimensionless variables and the London 
gauge div A = for the vector potential A , we write the 
system of GL equations in the following form 



-i^^n - ^1 * = * 



(l-|*l'), 



^sdA = -^5{z) J 



2" \^) J 2Di 



(1) 
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where 



J 2D 



= ^(**^2i)1'-*^2D*)-|*|'^, (3) 



is the density of superconducting current. The super- 
conducting wavefunction satisfies the following boundary 
conditions 
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and A — ^Hoplt^ far away from the superconductor. 
Here the distance is measured in units of the coherence 
length ^, the vector potential in ch/2e^, and the magnetic 
field in Hc2 — ch/2eS^^ — k^/2Hc. The ring is placed in 
the plane {x,y), the external magnetic field is directed 
along the z-axis, and the indices 2D, 3D refer to two- 
and three-dimensional operators, respectively. 

The giant vortex state in a circular configuration is 
characterized by the total angular momentum L through 
^ = ^J (p) ex.p{iL(j)), where p and (f> are the cylindrical 
coordinates. L is the winding number and gives the vor- 
ticity of the system. An arbitrary superconducting state 
is generally a mixture of different angular harmonics L 
due to the non-linearity of the GL equations. Never- 
theless, we can introduce an analog to the total angular 
momentum L which is still a good quantum number and 
which is in fact nothing else then the number of vortices 
in the system. 

To solve the system of Eqs.(|]-||), we generalized the 
approach of Ref. Q for disks to our fat ring configura- 
tion. We apply a finite-difference representation of the 
order parameter and the vector potential on a uniform 
Cartesian space grid (x,y), with typically 128 grid points 
over a distance 2Ro (i.e. the diameter of the ring), and 
use the link variable approach |]l9| , and an iteration pro- 
cedure based on the Gauss-Seidel technique to find ^. 
The vector potential is obtained with the fast Fourier 
transform technique where we set A\^^^fi^^\y\^ji^ = 
Hq {x, —y) /2 at the boundary of a larger space grid 
(typically Rs = ARo). 

To find the different vortex configurations, which in- 
clude the metastable states, we search for the steady- 
state solutions of Eqs.dBH) starting from different ran- 
domly generated initial conditions. Then we in- 
crease/decrease slowly the magnetic field and recalculate 
each time the exact vortex structure. We do this for each 
vortex configuration in a magnetic field range where the 
number of vortices stays the same. By comparing the 
dimensionless Gibbs free energies of the different vortex 
configurations 



F^y- 



2M- Ao • J2d-I*l 



(5) 



2i 



where integration is performed over the sample volume 
V, and Ag is the vector potential of the uniform mag- 
netic field, we find the ground state. The dimensionless 
magnetization, which is a direct measure of the expelled 
magnetic field from the sample, is defined as 



M = 



(g) - gp 
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(6) 



where Hq is the appUed magnetic field. {H) is the mag- 
netic field averaged over the sample or detector surface 
area S. 



III. SMALL RINGS: THE GIANT VORTEX 
STATE 

First we discuss a small superconducting ring. Al- 
though the system is circular symmetric in general, we 
are not allowed to assume that ^ {~p) — F [p) e*^^ be- 
cause of the nonlinear term in the GL-equations (|l|-||). 
Nevertheless, we found that for sufficiently small rings 
the confinement effects are dominant and this imposes 
a circular symmetry on the superconducting condensate. 
If we can assume axial symmetry, the dimensions of the 
GL equations will be reduced and thus the accuracy and 
the computation time will be improved. Therefore we 
will use this method for small rings. For a disk, this kind 
of solution with fixed L has been called 2D solution by 
Deo et al 0] in contrast to the 3D solution of the general 
problem. 

The ground state free energy i<" of a superconducting 
disk with radius Ro — 2.0^ and thickness d — 0.005^ and 
GL parameter k = 0.28 is shown in Fig. g for a hole 
in the center with radius Ri/S, = 0.0, 0.5, 1.0 and 1.5, 
respectively. The situation with R^ = corresponds to 
the situation of a superconducting disk without a hole 
which was already studied in Ref. H. With increasing 
hole radius Ri , we find that the superconducting/normal 
transition (this is the magnetic field where the free en- 
ergy equals zero) shifts appreciably to higher magnetic 
fields and more transitions between different i-states are 
possible before superconductivity disappears. Further- 
more, the thin dotted curve gives the free energy for a 
thicker ring with thickness d = 0.1^ and Ri = 1.0^. In 
comparison with the previous results for d = 0.005^, the 
free energy becomes more negative, but the transitions 
between the different L-states occur almost at the same 
magnetic fields. Thus increasing the thickness of the ring 
increases superconductivity which is a consequence of the 
smaller penetration of the magnetic field into the super- 
conductor. 

Experimentally, using magnetization measurements 
one can investigate the effect of the geometry and the 
size of the sample on the superconducting state. One 
of the main advantages of these measurements, in com- 
parison with resistivity measurements, is that they are 
non invasive and no conduction probes are needed which 
may distort the superconducting state. To investigate the 
magnetization of a single superconducting disk Geim et 
al [g^ used sub-micron Hall probes. The Hall cross acts 
like a magnetometer, where in the ballistic regime the 
Hall voltage is determined by the average magnetic field 



through the Hall cross region [£1|. Hence, by measur- 
ing the Hall resistance, one obtains the average magnetic 
field, and consequently, the magnetic field expelled from 
the Hall cross (Eq. (^)), which is a measure for the mag- 
netization of the superconductor m . As in the case of a 
superconducting disk 0] , the field distribution in the case 
of thin superconducting rings is extremely non-uniform 
inside as well as outside the sample and therefore the 
detector size will have an effect on the measured mag- 
netization. To understand this effect of the detector, we 
calculate the magnetization for a superconducting ring 
with outer radius Ro = 2.0^, thickness d = 0.005^ and 
two values of the inner radius Ri = 0.5^ and Ri — 1.5^ 
by averaging the magnetic field (see Eq. (^)) over several 
detector sizes S. The results are shown in Fig. ||(a,b) as 
a function of the applied magnetic field for Ri — 0.5^ and 
Ri — 1.5^, respectively. The solid curve (curve (i)) shows 
the calculated magnetization if we average the magnetic 
field over the superconducting volume, and curve (ii) is 
the magnetization after averaging the field over a circu- 
lar area with radius Ro, i.e. superconductor and hole. 
Notice that in the Meissner regime, i.e. Hq/Hc2 < 0.7 
where L = 0, the magnetization of superconductor + hole 
is larger than the one of the superconductor alone which 
is due to fiux expulsion from the hole. For L > 1 the 
reverse is true because now flux is trapped in the hole. 
Experimentally, one usually averages the magnetic field 
over a square Hall cross region. Therefore, we calculated 
the magnetization by averaging the magnetic field over 
a square region (curve(iii)) with width equal to 2i?o, i.e. 
the diameter of the ring. Curve (iv) shows the magne- 
tization if the sides of the square detector are equal to 
(2 + l/2)i?o. Increasing the size of the detector, decreases 
the observed magnetization because the magnetic field is 
averaged over a larger region which brings (H) closer to 
the applied field Hq 

Having the free energies of the different giant vortex 
configurations for several values of the hole radius varying 
from Ri ~ 0.0^ to i?i — 1.8^, we construct an equilibrium 
vortex phase diagram. Fig. H shows this phase diagram 
for a superconducting disk with radius Ro — 2.0^, thick- 
ness d — 0.005^ and for n = 0.28. The solid curves in- 
dicate where the ground state of the free energy changes 
from one L-state to another and the thick solid curve 
gives the superconducting/normal transition. The latter 
exhibits a small oscillatory behavior which is a conse- 
quence of the Little-Parks effect. Notice that the super- 
conducting/normal transition is moving to larger fields 
with increasing hole radius Ri and more and more flux 
can be trapped. In the limit Ri — > Ro, the critical mag- 
netic field is infinite and there are an infinite number 
of L-states possible which is a consequence of the en- 
hancement of surface conductivity for very small sam- 
ples y,|2^ . Because of the finite grid, we were not able 
to obtain accurate results for Ri « Ro- The dashed 
lines connect our results for hole radius Ri — 1.8^ with 
the results for Ri ^ Ro ||], where the transitions be- 
tween the different L-states occur when the enclosed flux 



is (j) = {L + l/2)(/)o, where 0o = ch/2e is the elemen- 
tary flux quantum. Notice that for rings of non zero 
width, i.e. Ri ^ Ro, the L ^ L + I transition occurs at 
higher magnetic field than predicted from the condition 
(j) = (L-|- 1/2)00- The discrepancy increases with increas- 
ing width of the ring and with increasing L. Starting 
from i?i = we find that with increasing Ri the Meiss- 
ner state disappears at smaller Hq. The hole in the center 
of the disk allows for a larger penetration of the magnetic 
field which favors the L = \ state. This is the reason why 
the L = ^ L = 1 transition moves to a lower exter- 
nal field while the L = \ ^ L = 2 transition initially 
occurs for larger Hq with increasing Ri. When the hole 
size becomes of the order of the width of one vortex the 
L = 1 ^r L ^ 2 transition starts to move to lower fields 
and the L — 2 state becomes more favorable. 

The effect of the thickness of the ring on the phase 
diagram is investigated in Fig. 0. The thin solid curves 
indicate where the ground state of the free energy changes 
from one L-state to another one, while the thick curve 
gives the superconducting/normal transition. Notice 
that the transition from the L = to the L = 1 state and 
the superconducting/normal transition depends weakly 
on the thickness of the ring. For increasing thickness d, 
the L = 2 state becomes less favorable and disappears for 
d > 0.7^. In this case, there is a transition from the L = 1 
state directly to the normal state. Increasing d stabilizes 
the different L-statcs up to larger magnetic fields. This 
is due to the increased expulsion of the applied field from 
the superconducting ring ^ . 

In the next step, we investigate three very important 
and mutually dependent quantities: the local magnetic 
field H, the Cooper pair density |^| and the current 
density j. We will discuss these quantities as a function 
of the radial position p. For this study we distinguish two 
situations, i.e. Ri <^ Ro and Ri < Ro- In the first case 
the sample behaves more like a superconducting disk, and 
in the second case like a superconducting loop. 

First, we consider a superconducting disk with radius 
Ro — 2.0^ and thickness d = 0.1^ with a hole with 
radius Ri — 0.5^ in the center. The free energy and 
the magnetization (after averaging over the superconduc- 
tor -|-hole) for such a ring are shown in Fig. |[ The dashed 
curves give the free energy and the magnetization for the 
different L-states, and the solid curve is the result for 
the ground state. Fig. |^(a,b) shows the local magnetic 
field H, Fig. |^(c,d) the Cooper pair density |*| , and 
Fig. |3(e,f) the current density j as function of the radial 
position p for such a ring at the L-states and magnetic 
fields as indicated by the open circles in Fig. ^(a,b). 

For low magnetic field, the system is in the L — state, 
i.e. the Meissner state, and the flux trapped in the hole 
is considerably suppressed. Hence, the local magnetic 
field inside the hole is lower than the external applied 
magnetic field as is shown in Fig. |7|(a) by curves 1 and 2. 
Please notice that the plotted magnetic field is scaled by 
the applied field Hq. In the L = state the superconduc- 



tor expels the magnetic field by inducing a supercurrent 
which tries to compensate the applied magnetic field in 
the superconductor and inside the hole. This is called the 
diamagnetic Meissner effect. As long as L equals zero, 
the induced current has only to compensate the mag- 
netic field at the outside of the ring and, therefore, the 
current flows in the whole superconducting material in 
the same direction and the size increases with increas- 
ing field. This is clearly shown in Fig. CTe) by curves 1 
and 2 where the current density j becomes more nega- 
tive for increasing Ho/Hc2- Notice also that the current 
density is more negative at the outside than at the inside 
of the superconducting ring which leads to a stronger de- 
pression of the Cooper pair density at the outer edge as 
compared to the inner edge of the ring (see curves 1 and 

2 in Fig. 0(c)). 

At Hq/Hc2 = 0.745, the ground state changes from 
the L = to the L — 1 state (see Fig. @(a)). Suddenly 
more flux becomes trapped in the hole (compare curve 1 
with curve 3 in Fig. ^a)), the local magnetic field inside 
the hole increases and becomes larger than the external 
magnetic field Hq. In the L = \ state, there is a sharp 
peak in the magnetic field at the inner boundary because 
of demagnetization effects. Consequently, more current 
is needed to compensate the magnetic field near the inner 
boundary than near the outer boundary (see curve 3 in 
Fig. n(e)). The sign of the current near the inner bound- 
ary becomes positive (the current direction reverses) , but 
the sign of the current near the outer boundary does not 
change. This can be explained as follows: Near the in- 
ner boundary [p > 1.0^) the magnetic field is compressed 
into the hole (paramagnetic effect), while near the outer 
boundary (p < 2.0^) the magnetic field is expelled to the 
insulating environment (diamagnetic effect). The sign re- 
versal of j occurs at p = p* and later we will show that 
the fiux through the circular area with radius p* is ex- 
actly quantized. At the L = to the L = 1 transition the 
maximum in the Cooper pair density (compare curves 1 
and 3 in Fig. 0(c)) shifts from p = Ri to p ^ Ro- Further 
increasing the external field increases the Cooper pair 
density near the inner boundary initially (compare curve 

3 and 4 in Fig. fflc)), because the flux in the hole has to 
be compressed less. The point p*, where j — 0, shifts 
towards the inner boundary of the ring. Further increas- 
ing the external magnetic fleld, the Cooper pair density 
starts to decrease (see curves 5 and 6 in Fig. ^(d)) and 
attains its maximum near the inner boundary. The cur- 
rent near the inner boundary becomes less positive (see 
curves 5 and 6 in Fig. 0(f)), i.e. less shielding of the 
external magnetic fleld inside the hole (see curves 5 and 
6 in Fig. 0(b)), and near the outer boundary j becomes 
less negative which shields the magnetic field from the 
superconductor-t-hole. Thus at the outer edge the local 
magnetic field has a local maximum which decreases with 
applied magnetic field Hq. 

At Hq/Hc2 ~ 2.0325, the ground state changes from 
the L = 1 state to the L = 2 state and extra flux is 
trapped in the hole. The changes in the magnetic fleld 



distribution, the Cooper pair density and the current 
density are analogous to the changes at the first transi- 
tion. For example, the magnetic field inside the hole in- 
creases compared to the external magnetic field (curve 7 
in Fig. n(b)), the radius p* increases substantially (curve 
7 in Fig. [7|(f)) and the maximum in the Cooper pair den- 
sity shifts to the outer boundary (curve 7 in Fig. 0(d)). 

For Ri <^ Ro and L > 0, the superconducting state 
consists of a combination of the paramagnetic and the 
diamagnetic Meissner state, like for a disk. For Ri < 
Ro we expect that the sample behaves like a loop and, 
hence, the superconducting state is a pure paramagnetic 
Meissner state or a pure diamagnetic Meissner state. 

We consider a superconducting disk with radius Ro — 
2.0^ and thickness d = O.lf with a hole with radius 
i?i = 1.8^ in the center. The free energy and the magne- 
tization (after averaging over the superconductor+hole) 
for such a ring are shown in Fig. H. The dashed curves 
give the free energy and the magnetization for the dif- 
ferent L-states, and the solid curve is the result for the 
ground state. Fig. §(a,b) shows the local magnetic field 
H, and Fig. p|(c,d) the current density j as function of the 
radial position p for such a ring at the L-states and mag- 
netic fields as indicated by the open circles in Fig. ||(a,b). 
The Cooper pair density has almost no structure and is 
practically constant over the ring and will, therefore, not 
be shown. 

For L = 0, the situation is the same as for Ri <C Ro- 
The magnetic field is expelled from the superconductor 
and the hole to the outside of the system, i.e. diamag- 
netic Meissner effect. The current flows in the whole 
superconducting material in the same direction (curves 
1 and 2 in Fig. 0(a)) and the size increases with in- 
creasing external field Hq (curves 1 and 2 in Fig. |^(c)). 
At Ho/Hc2 — 0.27, the ground state changes from the 
L = state to the L — I state and suddenly more flux is 
trapped in the hole. The local magnetic field inside the 
hole becomes larger than the external field Hq and there 
is a sharp peak near the inner boundary (curves 3 and 4 in 
Fig. 11(a)). In contrast to the situation for Ri <$^ Ro, there 
is no peak near the outer boundary, which means that the 
magnetic field is only expelled to the hole, i.e. paramag- 
netic Meissner effect. The induced current flows in the 
reverse direction in the whole superconductor (curves 3 
and 4 in Fig. ^(c)). For increasing external magnetic 
field, the magnetic field inside the hole, the height of the 
demagnetization peak and hence the size of the current 
decrease (see curves 3 and 4 in Fig. y(a,c)). Further in- 
creasing the field, the superconducting state transforms 
into a diamagnetic Meissner state. The magnetic field is 
now expelled to the outside of the sample (curves 5 and 6 
in Fig. ^(b)) and the direction of the current is the same 
everywhere in the ring (curves 5 and 6 in Fig. ||(d)). At 
Hq/Hc2 = 0.82, the ground state changes from the L — 1 
state to the L = 2 state. The changes in the magnetic 
field distribution (see curve 7 in Fig. ^(b)) and the cur- 
rent density (see curve 7 in Fig. ^(c)) are analogous to 
the changes at the first transition. The diamagnetic state 



transforms into a paramagnetic state. 

For a narrow ring with finite width, the superconduc- 
tor is in the paramagnetic or the diamagnetic Meissner 
state, like for a superconducting loop. Contrary to this 
infinitely narrow ring case, for narrow finite width rings 
the superconducting state can also consist of a combina- 
tion of these two states, i.e. the direction of the super- 
current in the inner part of the ring is opposite to the 
outer part. 

Now we will investigate the flux quantization in the 
fat ring. Fig. nfl shows the flux through the hole (a) 
and through the superconductor-|-hole (b) as function of 
the applied magnetic field Hq for a superconducting disk 
with radius Ro = 2.0^ and thickness d = 0.1^ with a hole 
in the center with radius Ri — 1.0^. The dashed curves 
show the flux for the different L-states, the solid curve 
for the ground state and the thin solid line is the flux 
through the hole if the sample is in the normal state. It 
is apparent that the flux through the hole (or through 
the superconducting ring -I- hole) is not quantized (see 
Fig. |0|(a)). At Ho/Hc2 = 0.4575 suddenly more flux 
enters the hole and the ground state changes from the 
L = state to the L = 1 state. At this transition also 
the flux increase through the hole is not equal to one flux 
quantum (j)^. It is generally believed that the flux through 
a superconducting ring is quantized. But as was shown in 
Ref. pq] this is even no longer true for hollow cylinders 
when the penetration length is larger than the thickness 
of the cylinder wall. The present result is a generalization 
of this observation to mesoscopic ring structures. Note 
that for the case of Fig. nO the penetration length is A/^ = 
0.28 and the effective penetration length A/^ = 0.78 is 
comparable to the width of the ring {Ro — Ri) /£, = 1.0. 

For L > 0, the superconducting current equals zero at 
a certain 'effective' radial position p = p*. It is the flux 
through the circular area with radius p* which is quan- 
tized and not necessarily the flux through the hole of our 
disk. Following Ref. ||l^ we integrate the supcrcurrent 
j over a closed contour C lying entirely inside the su- 
perconducting material which embraces the opening, one 
finds 
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K^r j + A ] -di = L(j)o. 



(4) 



Starting from the definition of flux and using Stokes' the- 
orem, the flux can be written as 



= H ■d~S = / rot~l ■ d~^ = 



A-dl 



(5) 



If the contour C is chosen in such a way that the current 
i ^ js ^ on this contour, then the fiux through the 
surface area bounded by this contour is quantized 



^V 



L(f)o. 



(6) 



To demonstrate that this is indeed true we show in 
Fig. |l^(a) the current density as a function of the radial 



position p and in Fig. [Tl|(b) the flux through a circular 
area with radius p for a superconducting disk with ra- 
dius Ro — 2.0^ and thickness d = 0.1^ with a hole in 
the center with radius Ri = 1.0^ in the presence of an 
external magnetic field Ho/Hc2 = 1.6075 for the case of 
three different giant vortex states; i.e. L — 1,2, 3. In the 
L = 1 state, the current density equals zero at a distance 
P*/S, ~ 1-16 from the center and the flux through an 
area with this radius is exactly equal to one flux quan- 
tum (j)o. For L = 2 and L = 3, the current density j 
equals zero at p* /£, ~ 1.56 and 1.91, respectively, and the 
flux through the area with this radius p* is exactly equal 
to 2^0 and 300, respectively. We find that p* depends 
on the external applied magnetic field and on the value 
of L, contrary to the results of Arutunyan and Zharkov 
|0,n who found that the effective radius p* is approx- 
imately equal to the geometric mean square of the in- 
ner radius Ri and the outer radius Ro of the cylinder; 
i.e. p* = [RrRo)^^'^ , which for the case of Fig,^|l|(a) 
would give p* = 1.41^. The results of Refs. |l^,|ljl were 
obtained within the London limit. The dependence of 
p* as function of the applied magnetic field is shown in 
Fig. ^(a). The dashed curves give the p* of the dif- 
ferent L-states. For increasing magnetic field and fixed 
L, the value of p* decreases, i.e. the "critical" radius 
moves towards the inner boundary. The solid curve gives 
p* for the ground state. At the L — *■ L -I- 1 transitions, 
p* jumps over a considerable distance towards the out- 
side of the superconducting ring, the size of the jumps 
decreases with increasing L. The dotted lines in the fig- 
ure give the two boundaries of the superconducting ring: 
the outer boundary Ro — 2.0^ and the inner boundary 
Ri = 1.0^. Remark that in Fig. 12 there is no p* given 
for the L = state, because only the external magnetic 
field has to be compensated so that the current has the 
same sign everywhere inside the ring and there exists no 
p* . In Fig. |2|(b) we repeated this calculation for a hole 
with radius Ri = 1.5^ where superconductivity remains 
to higher magnetic fields and many more L ^r L+1 tran- 
sitions are possible. The result of Refs. P^Jl^l gives in 
this case p* = 1.73^. Remark that the results for the 
L — \ and the L — 2 states are not connected. The rea- 
son is that just before the L = 1 ^ L = 2 transition the 
critical current for the L = \ state is strictly positive in 
the whole ring and hence p* is not defined. Notice that 
the results in Fig. [l^(a,b) oscillate around the average 
value p* = ^RiRo as given by Refs. ]l3| , p^ . 

For even narrower rings the current in the ring is 
mostly diamagnetic or paramagnetic and the region in 
magnetic field over which both directions of current oc- 
cur in the ring becomes very narrow. This is illustrated 
in Fig. O for a ring with outer radius Ro — 2.0^ and inner 
radius Ri — 1.8f . We plot in Fig. 0(a) the flux through 
the hole which becomes very close to the external flux, 
i.e. the flux without any superconductor. In Fig. |l3|(b) 
the value of the current at the inner and the outer side of 
the ring is shown which illustrates nicely that over large 



ranges of the magnetic field the current in the ring flows 
in one direction. The free energy becomes minimum (see 
Fig. tt3(c)) at a magnetic field where the current in the 
rings flows in both directions. 



IV. LARGE RINGS: THE MULTI- VORTEX 
STATE 

Until now, we considered only small rings. In such 
rings, the confinement effect dominates and we found 
that only the giant vortex states are stable and possible 
multi- vortex states have always larger energies if they ex- 
ist. Now we will consider larger superconducting disks in 
which multi-vortex states can nucleate for certain mag- 
netic fields. As an example we take a fat ring with outer 
radius Ro — 4.0^, thickness d = 0.005^, k = 0.28 and for 
different values of the inner radius. Fig. tt4(a) shows the 
free energy for such a ring with inner radius Ri — 0.4^ as 
a function of the applied magnetic field Hq . The differ- 
ent giant vortex states are given by the thin solid curves 
and the multi-vortex state by the thick solid curves. The 
open circles indicate the transitions from the multi- vortex 
state to the giant vortex state. In this ring, multi-vortex 
states exist with winding number L = 3 up to 7. For 
i = 3, 4, 5 multi- vortices occur both as mctastable states 
as well as in the ground state, while for L ~ 6,7 they are 
only found in the mctastable state. Notice there is no dis- 
continuity in the free energy at the transitions from the 
multi- vortex state to the giant vortex state for fixed wind- 
ing number L. Fig. H(b) shows the magnetization M for 
this ring as a function of Hq after averaging the field H 
over the superconducting ring (without the hole). The 
dashed curves give the results for the different L-states, 
the thin solid curve for the ground state, the thick solid 
curves for the multi-vortex states and the open circles 
indicate the transition from the multi- vortex state to the 
giant vortex states. Notice that the latter transitions are 
smooth, there are no discontinuities in the magnetization. 

Now, we investigate the flux (j) through the hole for the 
L — 4 multi- vortex and giant vortex state for the case of 
the above ring. Fig. UB shows the flux through a cir- 
cular area of radius p for different values of the applied 
magnetic field Hq. Curves 1,2,3 and 4 are the results for 
Hq/Hc2 = 0.72, 0.795, 0.87 (i.e. multi-vortex states) and 
0.945 (i.e. giant vortex state), respectively. There is no 
qualitative difference between the 4 curves, i.e. no quali- 
tative difference between the multi- vortex states and the 
giant vortex state. In the inset, we show the flux through 
the hole with radius Ri — 0.4^ and through the super- 
conductor -|-hole as a function of the applied magnetic 
field for a fixed value of the winding number, i.e. L = 4. 
The solid circles indicate the magnetic fields considered in 
the main figure and the open circle indicates the position 
of the transition from multi-vortex state to giant vortex 
state. The slope of the curves increases slightly at the ap- 
plied magnetic field, where there is a transition form the 
multi-vortex state to the giant vortex state. This agrees 



with the result for a disk ||] that the giant ^->multi- vortex 
transition is a second-order phase transition. 

The Cooper pair density h/;| for the previous four con- 
figurations is shown in Fig. ^. The darker the region, the 
larger the density and thus vortices are given by white re- 
gions. At the magnetic field Ho/Hc2 = 0.72 (Fig. 0(a)) 
we see clearly three multi- vortices. With increasing mag- 
netic field, these multi- vortices start to overlap and move 
to the center (Fig. nQ(b,c)) and finally they combine to 
one giant vortex in the center (Fig. 0(d)) Please notice 
that a theory based on the London limit will not be able 
to give such a complicated behaviour, because in such a 
theory vortices are rather point-like objects. 

The free energies of the different vortex configurations 
were calculated for different values of the hole radius 
which we varied from i?^ = to Ri = 3.6^. From these 
results we constructed an equilibrium vortex phase dia- 
gram. First, we assumed axial symmetry, where only gi- 
ant vortex states occur and the order parameter is given 
by * (7?) = F (p) e'-^". In the phase diagram (Fig. |l^) 
the solid curves separate the regions with different num- 
ber of vortices (different L-states) . In the limit Ri ^ 
we find the previous results of Ref. p| for a superconduct- 
ing disk. The radius of the giant vortex Rl in the center 
of the disk increases with increasing L, because it has 
to accommodate more flux, i.e. Rl/S, ~ \JL/{Hq/Hc2)- 
Hence, if we make a little hole in the center of the disk, 
this will not influence the L ~^ L + 1 transitions as long 
as Ri <^ Rl as is apparent from Fig. Hj. For sufficient 
large hole radius Ri, the hole starts to influence the giant 
vortex configuration and the magnetic field needed to in- 
duce the L -^ L+1 transition increases. For example, the 
transition field from the L — 7 state to the L = 8 state 
reaches its maximum for a hole radius Ri « 2.0^ which 
occurs for Hq/Hc2 ~ 1-5. The above rough estimate gives 
Rl/£, ^ 2.16 which is very close to Ri/^ ~ 2.0. Further 
increasing Ri, the hole becomes so large that more and 
more flux is trapped inside the hole, and consequently 
a smaller field is needed to induce the L ^ L + \ tran- 
sition. Because of flnite grid size, we were limited to 
Ri < 3.6^. The results we flnd for Ri = 3.6^ are ex- 
trapolated to = {L + 1/2)00 for Ri = Ro- The thick 
curve in Fig. ^ gives the superconducting/normal tran- 
sition. For low values of Ri, this critical magnetic field 
is independent of Ri, because the hole is smaller than 
the giant vortex state in the center and hence the hole 
does not influence the superconducting properties near 
the superconducting/normal transition. For Ri > 2.0^ , 
this field starts to increase drastically. Therefore, more 
and more L-states appear. In the limit Ri — + Ro, the 
critical magnetic field is infinite and there are an infinite 
number of L-states possible which is a consequence of 
the enhancement of surface conductivity for very small 
superconducting samples |22|. 

Next, we consider the general situation where the or- 
der parameter is allowed to be a mixture of different giant 
vortex states and thus we no longer assume axial symme- 
try of the superconducting wavefunction. We found that 



the transitions between the different L-states are not in- 
fluenced by this generalization, but that for certain mag- 
netic fields the ground state is given by the multi-vortex 
state instead of the giant vortex state. In Fig. |l^ the 
shaded regions correspond to the multi- vortex states and 
the dashed curves are the boundaries between the multi- 
vortex and the giant vortex states. For L = 1, the single 
vortex state and the giant vortex state are identical. In 
the limit Ri — > 0, the previous results of Ref. Q for a 
superconducting disk are recovered. For increasing hole 
radius Ri, the L = 2 multi- vortex state disappears as a 
ground state for Ri > 0.15^. If Ri is further increased, 
the ground state for L = 5 up to L = 9 changes from 
giant vortex state to multi-vortex state and again to gi- 
ant vortex state. For example, for Ri = 2.0^ the multi- 
vortex state exists only in the L = 9 state. Notice that 
for small Ri the region of multi-vortex states increases 
and consequently the hole in the center of the disk sta- 
bilizes the multi-vortex states, at least for L > 2. For 
large Ri, i.e. narrow rings, the giant vortex state is the 
energetic favorable one because confinement effects start 
to dominate which impose the circular symmetry on tp. 
For fixed hole radius Ri < 2.0^ and increasing magnetic 
field we find always at least one transition from giant 
vortex state to multi-vortex state and back to giant vor- 
tex state (re-entrant behaviour). Remark that the ground 
state for L > 10 is in the giant vortex state irrespective 
of the value of the magnetic field. Near the supercon- 
ducting/normal transition the superconducting ring is in 
the giant vortex state because now superconductivity ex- 
ists only near the edge of the sample and consequently 
the superconducting state will have the same symmetry 
as the outer edge of the ring and thus it will be circular 
symmetric. 

Finally, for L > 3 the multi-vortex state not necessar- 
ily consists of L vortices in the superconducting mate- 
rial. Often they consist of a combination of a big vortex 
trapped in the hole and some multi-vortices in the su- 
perconducting material. This is clearly shown in Fig. |l8| 
where a contour plot of the local magnetic field is given 
for a superconducting disk with radius Ro — 4.0^ and 
thickness d = 0.005^ with a hole in the center with ra- 
dius i?, = 0.6^ (Fig. |l|(a)) and R, = 1.0^ (Fig. |l|(b)). 
As usual we took k = 0.28. The dashed thick circle 
is the outer radius, while the small solid black circle is 
the inner radius of the ring. Low magnetic fields are 
given by light regions and darker regions indicate higher 
magnetic fields. In this way, multi-vortices in the su- 
perconducting area are dark spots. 



In Fig. |l8|(a) the lo- 
cal magnetic field is shown for an applied magnetic field 
Ho/Hc2 — 0.895. Although the winding number is L = 4, 
there are only 3 vortices in the superconducting material 
and one vortex appears in the hole in the center of the 
disk. This is clearly shown in Fig. |9|(a) where the phase 
if of the order parameter is shown along different cir- 
cular loops T* — > Ce^^ inside the superconductor. The 
solid curve gives the phase near the outer edge of the 
ring (C = 3.95^) and the dashed curve near the inner 



edge of the ring (C = 0.7^). When encirchng the ring, 
the phase difference A(^ in the first case is 4 times 27r, 
while in the second case it is /S.(p = f x 2tt. The phase 
difference is always given by Atp — L x 27r, with L the 
winding number. In Fig. f 8[b) a contour plot of the local 
magnetic field is shown for a ring with _R, = 1.0^ which 
leads to L = 6 at Hq/Hc2 — l.f45. Only 4 vortices are 
in the superconducting material and one giant vortex in 
the center (partially in the hole) with L = 2. Fig. p^(b) 
shows the phase (/? of the order parameter for C — 3.95^ 
(solid curve) and 1.2^ (dashed curve), where the phase 
differences are At^ = 6 x 27r and 2 x 27r, respectively. No- 
tice that the flux through the hole equals </> ~ O.190o 
for the case of Fig. [l8|(a) and </> ~ O.360o for the case of 
Fig. |fj(b) and is thus not equal to a multiple of the flux 
quantum 0o- 



V. NON-SYMMETRIC GEOMETRY 

So far, we investigated the influence of the size of 
the hole on the vortex configuration for superconducting 
rings of different sizes. We found that for small rings, 
only the axially symmetrical situation occurs, i.e. the gi- 
ant vortex states. For large rings the multi-vortex state 
can be stabilized for certain values of the magnetic field. 
In the next step, we purposely break the axial symmetry 
by moving the hole away from the center of the super- 
conducting disk over a distance a. 

As an example, we consider a superconducting disk 
with radius Ro = 2.0^ and thickness d = 0.005^ with 
a hole with radius Ri = 0.5^ moved over a distance 
a = 0.6^ in de x-direction. Fig. EO^ shows the free en- 
ergy and magnetization (defined through the field ex- 
pelled from the superconducting ring without the hole) 
as function of the magnetic field. The solid curve in- 
dicates the ground state, while the dashed curves indi- 
cate the metastable states for increasing and decreas- 
ing field. The vertically dotted lines separate the re- 
gions with different winding number L. Notice that 
hysteresis is only found for the first transition from the 
i = to the L = 1 state and not for the higher tran- 
sitions which arc continuous. At the transition from 
the L = 1 state to the L — 2 state and further to the 
L — 3 state, the free energy and the magnetization vary 
smoothly (see Fig. ^(a,b)). In the insets of Fig. pO|(a), 
we show the Cooper pair density for such a sample at 
magnetic field Hq/Hc2 = 0.145, 1.02, 2.145 and 2.52, re- 
spectively, where the ground state is given by a state 
with L — 0,1,2,3 respectively. High Cooper pair den- 
sity is given by dark regions, while light regions indicate 
low Cooper pair density. For Ho/Hc2 — 0.145, we find 
a high Cooper pair density in the entire superconduct- 
ing ring. There is almost no flux trapped in the circular 
area with radius smaller than Rq. After the first tran- 
sition at Ho/Hc2 ~ 0.75, suddenly more flux is trapped 
in the hole which substantially lowers the Cooper pair 



density in the superconductor. Notice that the trapped 
flux tries to restore the circular symmetry in the Cooper 
pair density and that the density of the superconducting 
condensate is largest in the narrowest region of the super- 
conductor. The next inset shows the Cooper pair density 
of the L = 2 state where an additional vortex appears. 
Some fiux is passing through the hole (i.e. winding num- 
ber is one around the hole) , while one flux line is passing 
through the superconducting ring and a local vortex (the 
normal region with zero Cooper pair density) is created. 
In the L — 3 state superconductivity is destroyed in part 
of the sample which contains flux with winding number 
L = 2 and the rest of the flux passes through the hole. 
Hence, by breaking the circular symmetry of the system, 
multi-vortex states are stabilized. Remember that for the 
corresponding symmetric system, i.e. a = 0, only giant 
vortex states were found. 

Having the magnetic flelds for the different L ^' L + 1 
transitions for superconducting rings with different posi- 
tions of the hole, i.e. different values of a, we constructed 
the phase diagram shown in Fig. |2l|. The thin curves 
(solid curves when the magnetization is discontinuous 
and dashed curves when the magnetization is continu- 
ous) indicate the magnetic fleld at which the transition 
from the L-state to the i-|-l-state occurs, while the thick 
solid curve gives the superconducting/normal transition. 

In order to show that the stabilization of the multi- 
vortex state due to an off-center hole is not peculiar 
for Ro = 2.0^, we repeated the previous calculation for 
a larger superconducting disk with radius Ro = 5.0^ 
and thickness d = 0.005^ containing a hole with radius 
Ri = 2.0^. In Fig. ^(a) the Cooper pair density is shown 
for such a system with the hole in the center, while in 
Fig. P3(b) the hole is moved away form the center over 
a distance a — 1.0^ in the negative y-direction. The ex- 
ternally applied magnetic field is the same in both cases, 
Ha/Hc2 ~ 0.77, which gives a winding number L = 4 
and L ~ 5 for the ground state of the symmetric and 
the non-symmetric geometry, respectively. The assign- 
ment of the winding number can be easily checked from 
Figs. g2|(c,d) which show contour plots of the correspond- 
ing phase of the superconducting wavefunction. If the 
hole is at the center of the disk, the ground state is a 
giant vortex state. Moving the center of the hole to the 
position {x/^,y/S,) = (0, —1), two vortices appear in the 
superconducting material while the hole contains three 
vortices. Notice that in this case, although the magnetic 
field is kept the same and the amount of superconduct- 
ing material is not altered, changing the symmetry of the 
system alters the winding number. 



VI. CONCLUSION 

In conclusion, we studied the superconducting state of 
thin superconducting disks with a hole. The effect of 
the size and the position of the hole on the vortex con- 



figuration was investigated. For small superconducting 
disks with a hole in the center, only giant vortex states 
exist and for increasing hole radius Ri more and more 
L-states occur before the superconductor becomes nor- 
mal. For larger superconducting disks with a hole in the 
center, we found multi-vortex states in a certain mag- 
netic field range. For certain fixed hole radius, and for 
increasing magnetic field, the giant vortex state changes 
into a multi-vortex state and back into the giant-vortex 
state (re-entrant behaviour) before superconductivity is 
destroyed. Near the superconducting/normal transition 
and for a narrow superconducting ring (i.e. Ri « Ro) we 
always found the giant vortex state as the ground state 
irrespective of the size, thickness and width of the ring. 
The effect of the position of the hole, i.e. decreasing the 
symmetry of the system, was also investigated. Moving 
the hole off-center: 1) can transform the L ^ L + 1 tran- 
sition into a continuous one, 2) the stability of metastable 
states is strongly reduced, 3) it favours the multi-vortex 
state even for small disks, and 4) the winding number L 
can increase even at a fixed magnetic field. 

The flux through the hole is not quantized. We were 
able to define an effective ring size p* such that inside 
this ring the flux is exactly quantized. The value of 
Ri < P* < Ro depends on L and is an oscillating function 
of the magnetic field. For narrow rings it is only possi- 
ble to define such a p* in narrow ranges of the magnetic 
field and the flux through the hole is very close to the 
applied flux. The magnetic fields from the screening cur- 
rents are too small to substantially modify the flux inside 
the ring. On the other hand, the magnetic field incre- 
ment AH or the flux increase Acp = ttR^AH to induce 
the L —^ L + I transition is only quantized for narrow 
rings. This is illustrated in Fig. g^ in case of Ro = 4.0^ 
for different values of the inner radius. Notice that for 
Ri ^ Ro we find that A^ is an oscillating function of L. 
It is substantially larger than 0o for small L, it is smaller 
than 00 for intermediate L and it approaches ^o from 
above for large L. We found earlier that for Ri w Ro 
the flux through p* = y/RiRo is quantized in 0o and 
therefore we expect A0* = n (p*)" AH = 0o- Our defini- 
tion of A(j) considers the flux through the superconduct- 
ing ring -|- hole which for the condition A0* = 0o gives 
A0/0O w 6.67,2.5, 1.54 and 1.11 for Ri/£_ = 0.6,1.6,2.6 
and 3.6, respectively. These results for Ri/£, = 2.6 and 
3.6 agree rather well with our numerical results presented 
in Fig. g^; i.e. A4>/(f>o « 1.5 and 1.1, respectively. The 
results presented in Fig. ^ agree qualitatively with the 
recent theoretical results of Bruyndoncx et al |1^ who 
studied rings in a homogeneous magnetic field. 
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FIG. 1. The configuration: a superconducting disk with 
radius Ro and thickness d with a hole inside with radius Ri 
which is placed a distance a away from the center. 



FIG. 2. The ground state free energy as a function of the 
applied magnetic field Ho of a superconducting disk with ra- 
dius Ro = 2.0^, thickness d = 0.005^ and k — 0.28 for a hole 
in the center with radius Rt/^ = 0.0, 0.5, 1.0 and 1.5 respec- 
tively. The thin dotted curve gives the free energy of a thicker 
ring with thickness d = 0.1^ and Ri = 1.0^. The free energy 
is in units of Fq = H'^V/8n. 



FIG. 3. The magnetization as a function of the applied 
magnetic field Hq for a superconducting disk with radius 
Ro = 2.0^ and thickness d = 0.005^ and k = 0.28 for a 
hole in the center with radius Ri/^ = 0.5 (a) and Ri/^ — 1.5 
(b). Curve (i) is the calculated magnetization if we average 
the magnetic field over the superconducting volume; curve 
(ii) after averaging the field over the area with radius Ro, i.e. 
superconductor and hole; curves (iii) and (iv) after averaging 
the magnetic field over a square region with widths equal to 
2Ro and (2 + l/2)i?o, respectively. 



FIG. 4. Phase diagram: the relation between the hole ra- 
dius Ri and the magnetic fields Hq at which giant vortex tran- 
sitions L ~* L + 1 takes place for a superconducting disk with 
radius Ro = 2.0^ and thickness d = 0.005^ and for k = 0.28. 
The solid curves indicate where the ground state of the free 
energy changes from one L-state to another one and the thick 
solid curve gives the superconducting/normal transition. 



FIG. 8. The free energy (a) and the magnetization after 
averaging over the superconductor-f hole (b) as a function of 
the applied magnetic field for a superconducting disk with 
outer radius Ro = 2.0^ and thickness d = 0.1^ with a hole 
with radius Ri = 1.8^ in the center (k — 0.28) for the differ- 
ent I/-states (dashed curves) and for the ground state (solid 
curves). The open circles are at: (1) L = 0, Ho/Hc2 = 0.27; 
(2) L = 0, Ho/Ho2 = 0.395; (3) L = 1, Ho/H,2 = 0.27; (4) 
L = 1, Ho/Ho2 = 0.395; (5) L = 1, Ho/H,2 = 0.695; (6) 
L = 1, Ho/Ho2 = 0.82; and (7) L = 2, H0/H02 = 0.82. 



FIG. 9. The local magnetic field H (a,b) and the current 
density 7 (c,d) for the situations indicated by the open circles 
in Fig. W as a function of the radial position p for a supercon- 
ducting disk with radius Ro — 2.0^ and thickness d = 0.1^ 
with a hole with radius Ri = 1.8^ in the center (k = 0.28). 



FIG. 10. The fiux through (a) the hole and (b) the super- 
conductor -I- the hole as a function of the applied magnetic 
field Ho for a superconducting disk with radius Ro = 2.0^ 
and thickness d = 0.1^ with a hole in the center with radius 
Ri = 1.0^. The dashed curves show the fiux for the different 
I/-states, the solid curve for the ground state. The thin solid 
line is the fiux in the absence of superconductivity. 



FIG. 11. The current density (a) and the fiux (b) as a func- 
tion of the radial position p for a superconducting ring with 
Ro = 2.0c, Ri = l.OC, d = O.IC and k = 0.28 for L = 1 (solid 
curves), L — 2 (dashed curves) and i = 3 (dotted curves) at 
an applied magnetic field Ho/Hc2 = 1.6075. 



FIG. 5. The same as Fig. ^ but now for varying thickness 
d of the ring for fixed Ro = 2.0^, Ri = 0.5$ and k = 0.28. 



FIG. 6. The free energy (a) and the magnetization after 
averaging over the superconductor-|-hole (b) as a function of 
the applied magnetic field for a superconducting disk with 
outer radius Ro = 2.0^ and thickness d = 0.1^ with a hole 
with radius Ri = 0.5^ in the center (k = 0.28) for the differ- 
ent L-states (dashed curves) and for the ground state (solid 
curves). The open circles are at: (1) L — 0, H0/H02 ~ 0.745; 
(2) L = 0, H0/H02 = 0.995; (3) L = 1, H0/H02 = 0.745; (4) 
L = 1, H0/H02 = 0.995; (5) L = 1, H0/H02 = 1.7825; (6) 
L = 1, H0/H02 = 2.0325; and (7) L = 2, H0/H02 = 2.0325. 



FIG. 12. The dependence of the effective radius p* as func- 
tion of the applied magnetic field for a superconducting disk 
with radius Ro = 2.0^ and thickness d = 0.1^ with a hole in 
the center with radius Ri = 1.0^ (a), and Ri — 1.5^ (b). The 
dashed curves show p* for the different L-states and the solid 
curve is for the ground state. 



FIG. 13. (a) The flux through the hole of the ring as func- 
tion of the applied magnetic field. The thin dashed line is the 
applied fiux. (b) The current at the inner side (solid curve) 
and at the outer side (dashed curve) of the ring, and (c) the 
free energy of the ground state as function of the applied 
magnetic field. 



FIG. 7. The local magnetic field H (a,b), the Cooper pair 
density I*!*! (c,d) and the current density j(e,i) for the situ- 



ations indicated by the open circles in Fig. ra as a function of 
the radial position p for a superconducting disk with radius 
Ro = 2.0^ and thickness d — 0.1^ with a hole with radius 
Ri = 0.5^ in the center (k = 0.28). 
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FIG. 14. (a) The free energy as function of the apphed 
magnetic field for a superconducting ring with Ro — 2.0^, 
Ri = 0.4^, d = 0.005^, and k = 0.28. The different gi- 
ant vortex states are shown by the thin solid curves and the 
multi-vortex states by the thick solid curves, (b) The mag- 
netization for the same sample after averaging over the su- 
perconducting ring only. The different L-states are given by 
the dashed curves, the ground state by the thin solid curves 
and the multi-vortex states by the thick solid curves. The 
transitions from multi-vortex state to giant vortex state are 
indicated by the open circles. 



FIG. 15. The flux (f) through a circular area of radius p for 
different values of the applied magnetic field; Hq. Curves 1,2,3 
and 4 give the results for Ho/Hc2 = 0.72 (1), 0.795 (2), 0.87 
(3) (i.e. multi-vortex states) and 0.945 (4) (i.e. giant vortex 
state), respectively. The inset shows the flux through the 
hole with radius Ri = 0.4^ and through the superconductor 
-I- hole as a function of the applied magnetic field for a fixed 
value of the winding number, i.e. L = 4. The solid circles 
indicate the magnetic fields considered in the main figure and 
the open circles indicate the transition from multi-vortex to 
giant vortex state. 



FIG. 16. The Cooper pair density corresponding to the four 
situations of Fig. ^ Ho/H^2 = 0.72 (a); Ho/Hc2 = 0.795 (b); 
Ho/Hc2 = 0.87 (c); and Hq/Hc2 = 0.945 (d). Dark regions 
indicate high density, light regions low density. The thick 
circles indicate the inner and the outer edge of the ring. 



FIG. 19. The phase tp of the order parameter calculated on 
a circle r -^ Ce^^ as a function of the angle x for a super- 
conducting ring with Ro = 4.0^, d = 0.005^ and n = 0.28; 
(a) Ri = 0.6C, Ho/H^2 = 0.895, L = 4 and C = 3.95^ 
(solid curve) and C = 0.7^ (dashed curve); (b) Ri = 1.0^, 
Ho/Hc2 = 1.145, L = 6 and C = 3.95^ (solid curve) and 
C = 1.2^ (dashed curve). 



FIG. 20. The free energy (a) and magnetization after av- 
eraging over the superconducting ring only (b) as function 
of the applied magnetic field for a superconducting disk with 
radius Ro — 2.0^ and thickness d — 0.005^ with a hole with 
radius Ri = 0.5^ moved over a distance a = 0.6^ in the 
a;-direction. The solid curve indicates the ground state, the 
dashed curve the results for increasing and decreasing field. 
The vertically dotted lines separate the regions with different 
vorticity L. The insets show the Cooper pair density at mag- 
netic field Ho/Hc2 = 0.145, 1.02, 2.145 and 2.52, where the 
ground state is given by a state with L = 0, 1, 2, 3 respectively. 
High Cooper pair density is given by dark regions, while light 
regions indicate low Cooper pair density. 



FIG. 21. The relation between the displacement a and the 
transition magnetic fields for a superconducting ring with 
i?o/C = 2.0, d/C = 0.005, Ri/^ = 0.5 and k = 0.28. The 
thin solid curves give the L —f L + 1 transitions when it 
corresponds to a first-order phase transition while a dashed 
curve is used when it is continuous. The thick solid curve is 
the superconducting/normal transition. 



FIG. 17. Equilibrium vortex phase diagram for a super- 
conducting disk with radius Ro — 4.0^, thickness d — 0.005^ 
with a hole with radius Ri in the center. The solid curves 
show the transitions between the different L-states, the thick 
solid curve shows the superconducting/normal transition and 
the dotted lines connect the results for Ri — 3.6^ with the 
results in the limit Ri — > Ro- The shaded regions indicate 
the multi-vortex states and the dashed curves separate the 
multi- vortex states from the giant vortex states. 



FIG. 22. The Cooper pair density for a superconducting 
disk with radius Ro = 5.0^ and thickness d = 0.005^ with a 
hole with radius Ri — 2.0^ (a) in the center, and (b) moved 
away form the center over a distance a = 1.0^ in the negative 
y-direction. The applied magnetic field is the same in both 
cases: Ho/Hc2 ~ 0.77. High Cooper pair density is given by 
dark regions, low Cooper pair density by light regions. The 
corresponding contour plot of the phase of the superconduct- 
ing wavefunction is given in (c) and (d). 



FIG. 18. Contour plot of the local magnetic field for a 
superconducting disk with radius Ro — 4.0^ and thickness 
d = 0.005^ (k = 0.28) with a hole in the center with (a) ra- 
dius Ri = 0.6C at Ho/Hc2 = 0.895 for L = 4; and (b) radius 
Ri = l.OC at Ho/Hc2 = 1.145 for L = 6. The dashed thick 
circle is the outer radius, the small solid thick circle the inner 
radius. Low magnetic fields are given by light regions and 
dark regions indicate higher magnetic fields. 



FIG. 23. The fiux increase A(/!> — tvRoAH needed to induce 
the L -^ L + 1 transition for a superconducting ring with 
Ro = 4.0^ for different values of the inner radius Ri. The 
interconnecting lines are a guide to the eye. 
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